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Abstract 

The appearance of a new phase of QCD, Quarkyonic Matter in the limit of large 
number of colors is studied within Nambu-Jona-Lassinio effective chiral model cou- 
pled to the Polyakov loop. The interplay of this novel QCD phase with chiral sym- 
metry restoration and color deconfinement is discussed. We find that at vanishing 
temperature and at large N c , the quarkyonic transition occurs at densities only 
slightly lower than that expected for the chiral transition. This property is also 
shown to be valid at finite temperature if the temperature is less than that of de- 
confinement. The position and iV c -dependence of chiral critical end point is also 
discussed. 
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1 Introduction 



The conventional view on the phase diagram of QCD is that high density 
strongly interacting matter is divided into two phases: the confined and the de- 
confined [TJ2 3 4 5 6 l71l8ll9lfT0|ll|12|fT3] . The phase diagram of QCD as a func- 
tion of baryon number chemical potential fi and temperature T as shown in 
Fig. 1, was originally envisioned by Cabibbo and Parisi [3] and has changed lit- 
tle conceptually since. Until very recently, the possible new physics in the QCD 
phase diagram were Color Superconducting phases which might be important 
at asymptotically high baryon number density and low temperatures [13]. In 
addition, a combination of efforts of Lattice Gauge Theory (LGT) [15] and 
effective model calculations |16II17|18II19II20|21] have given new insight on the 
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Fig. 1. The "phase diagram" of QCD presented as envisaged by Cabibbo and Parisi. 

position, the order and the universal properties [22] of the QCD phase dia- 
gram. 

It has been argued recently that there may be an additional phase, the Quarky- 
onic Phase, of dense QCD [23] 0- This phase was rigorously shown to exist in 
the limit of a large number of colors N c . In this limit, both the exponential of 
the free energy of a heavy test quark added to the system 

e-™ = ±(L), (1) 

and the baryon number density are order parameters [23] . The baryon number 
is an order parameter since (Nb) ~ e -/3( A ^s-Ms)_ Thus, for temperatures of 
T ~ Aqcd and with M B ~ iV c the (N B ) ~ e~ KNc -> at large N c , as long as 
the baryon number chemical potential is small compared to the baryon mass, 
i.e. \i <C Mb- When fi B > Mg, then baryons begin to populate the system 
and the baryon number density is non-zero. In Ref. [23 J, it was argued that 
there are at least three phases in QCD at large N c : the mesonic-phase which 
is confined and has zero baryon number density, the de-confined phase which 
has finite baryon number density, and the quarkyonic-phase which has finite 
baryon number density and is confined. The role of the chiral phase transition 
was not established. 

The reason for the existence of the quarkyonic world was because for any finite 
value of chemical potential for quarks, \iq = fi B /N c , quark loops do not affect 
the confining potential. The de-confinement temperature is at some T c and is 
independent of \iq. Therefore when baryons are added to the system, one can 
compress the baryons to very high chemical potential compared to Aqcd, and 
the baryonic matter remains confined. When \xq > \/W c Kqcdi then the effects 
of the quark loops are felt on the potential and there is de- confinement, but 

1 Speculations about related phases of matter were made in early strong coupling 
lattice studies of QCD at high baryon number density |24l25j . 
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Fig. 2. The phase diagram of QCD in large N c . We do not display either the Chiral 
or Color Superconducting phases on this plot. 
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Fig. 3. A hypothetical phase diagram including 1/N C effects (again ignoring the 
effects of the Chiral transition and Color Superconductivity). 

as A^c ^ oo, the density at which this occurs approaches infinity. 



The mesonic world is confined and has an energy density which scales as 
0(1) in powers of N c . The de-confined energy density scales as N%, due to 
unconfined gluons. The energy density of the quarkyonic world scales as N c , 
since both for baryonic matter and quark matter, the energy density is of order 
N c . The quarkyonic world may be visualized as a quasi free degenerate Fermi 
gas of quarks in a sea of thermally excited mesons and glueballs. The effects of 
confinement are important for quark interactions only near the Fermi surface. 
The bulk interactions deep inside the Fermi sea, even though in a confined 
phase, are described by perturbation theory. The name quarkyonic was chosen 
since it is a combination of baryonic and quark matter, and expresses the Yin- 
Yang nature of the matter. A hypothetical phase diagram of QCD in the large 
N c limit is shown in Fig. [2J 

In this paper we outline a theory which allows for explicit computation in the 
context of the PN JL model of QCD [TSlfTU] . This provides a concrete descrip- 
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Fig. 4. A guess for the phase diagram of QCD for realistic value of N c (without the 
Chiral and Color Superconducting phase transitions). 



tion of quarkyonic matter along with spontaneous chiral symmetry breaking 
and the effects of finite N c on the quarkyonic phase transition. In the large N c 
limit, we consider an exactly solvable model, which has the features expected 
for quarkyonic matter. We will argue that as the baryon number density in- 
creases from zero, the confinement- deconfinement phase transition weakens, 
in accord with the arguments of de Forcrand and Philipsen [26] . We also argue 
as one increases the temperature from zero, the quarkyonic phase transition 
weakens. The chiral phase transition, which is very close to that of the quarky- 
onic phase transition is first order, and is situated almost atop the region where 
the quarkyonic phase transition took place, until the critical temperature is 
reached where it then follows the deconfinement phase transition. 

As we decrease N c from asymptotically large values, there is some point where 
the first order deconfinement transition weakens and begins to disappear, as 
shown in Fig. 3. Eventually the low density confinement-deconfinement part of 
the phase transition completely disappears, leaving a line of first order phase 
transitions and a critical point which is a remnant of the chiral-quarkyonic 
phase transition, as shown in Fig. 4. These distinct branches of the phase 
diagram were shown previously in the work within the PNJL model [27f28] . 
and by Miura and Ohnishi in strong coupling lattice gauge theory |29j. Strictly 
speaking, the quarkyonic phase transition becomes a sharp cross over at finite 
N c . We find that in large N c , the quarkyonic transition occurs at densities 
slightly lower than that of the chiral transition. This difference in density is 
however so small that it may be an artifact of the model. As shown on Fig. 4, 
the cross over associated with the deconfinement phase transition continues 
on to higher densities, and becomes a separate line. This will be explained in 
detail later. 
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2 The PNJL Model in Large N c 



In order to study the QCD phase diagram in large N c we construct a chiral 
model where constituent quarks [5U] couple to effective gluon degrees of free- 
dom. Here we follow [TBfl8lll9| to introduce an extended Nambu-Jona-Lasinio 
model with Polyakov loops (PNJL model) 0. 

We take the Lagrangian for a constituent quark field ip as 

£ = ^(z 7 ^ M ~m + ifijo) V> + 1 {(W 2 + (^r 75 ^) 2 } -U($ [A], $L4]) , (2) 

where m is the current quark mass, /i is the quark chemical potential and f 
are Pauli matrices. In the following we restrict our discussion to two quark 
flavors, Nf = 2. An extension of the model to Nf > 2 is straightforward. 

The interaction between the quarks and the effective gluon field is implemented 
through a covariant derivative 

D ll = d li -iA ll , — d^A , (3) 

where A^ = gA a ^. Here g is the color SU(3) gauge coupling constant and A a 
are the Gell-Mann matrices. In the PNJL model, the transverse components 
of A are integrated out. If we assume that these reflect short distance degrees 
of freedom, the effective potential in terms of $ and $ should be at most 6'th 
order in the fields ($, $), since these are the relevant operators for the three 
dimensional space on which $ and $ exist 0. The variable $ is defined as a 
trace of the Wilson line L = Pexp |z f$ drAo(x, r)| in color space: 

$ = i-TrL , (4) 



and $ is the complex conjugate. 

We will work in mean field approximation for the fields $. The potential U 
must respect the Z(N C ) symmetry. For large N c , the most general potential is 
of the form [16J, 

^** + !(») 2 +^(**) 3 ), (5) 




One should keep in mind that the PNJL model describes statistical suppression 
of colored one- and two-quarks contributions which imitates color confinement. 
3 When we have a kinetic energy term (<9<I>) 2 of a scalar field we would have 
to rescale $ in such a way that it has the correct canonical dimension (d — 2)/2. 
Thus, the highest order term with a renormalizable coupling involves the 6th order 
in scalar fields. 
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a2 


6 4 


6 6 C 


0.787 


0.333 


-1.13 


-0.213 


1.00 5.35 



Table 1 

Set of parameters for the Polyakov-loop effective potential. 

with an overall constant C. One factorizes the N^ — l dependence to get at high 
T the pressure of an ideal gluon gas. The coefficients b{ must be chosen so that 
at high temperatures, there is spontaneous breaking of the Z(N C ) symmetry, 
and at low T the symmetry is restored. For N c > 3, the pure gauge theory 
has a first order phase transition corresponding to the de-confinement at some 
temperature T . At T = T one finds 

HTo) = -§ e . (6) 

We assume that 62 (T) has the following temperature dependence: 

T \ /TV 2 



b 2 (T) = a + 0l J + a 2 j , (7) 

with constant parameters <2j. The Polyakov loop expectation value ($) must 
be unity at asymptotically high temperature T — > 00. This leads to 

a = 64 + b 6 . (8) 

We fix the parameters Oj, fej and C taking 6 6 = 1 in such a way that (JSJ) 
describes the LGT observations for $(T ) and pressure P{T) = —U(T) in 
SU(N C = 3) pure gauge theory [31]. We list the resulting parameters in Ta- 
ble [H We assume that T is approximately independent of N c in U. This is 
supported by a finding in lattice QCD for several N c [32] where the decon- 
finement transition temperature approaching from low temperature phase is 
parameterized as 

T (N C ; n = 0) = ^0.596 + , (9) 

with the string tension a being independent of N c . Thus at large N c , T Q is 
Oil). 

In this paper we are interested in structural features of the PNJL model. Here 
we simply assume that there is a first order transition in the absence of quarks, 
and then determine the effect of quarks on this transition. Note that in the 
large N c limit, to leading order, the quarks do not affect the effective potential 
U. We will therefore assume a rigid background of fields $ and explore the 
consequences for fermion-induced phase transitions. 

Note also, that we need to specify which root of $ we take when the Z(N C ) 
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symmetry is broken. We take the real positive root, the root which is selected 
if there is a small breaking of the Z(N C ) symmetry induced by fermions. 

Using bosonization, this Lagrangian can be re-expressed as 

2 2 

C= - U- ° -zTrlnff- 1 . (10) 

2G 

with quark propagator 

S- 1 =i 7 ^- 7o A°-Af, (11) 
and dynamical quark mass 

M = m - (cr + ry 5 f • 7?) . (12) 
The field a in the mean field approximation, 

(a) = G(^) , (13) 
gets a nonzero expectation value from solving the gap equation. 

3 The NJL sector 



In the following we first consider the solution of the fermionic sector of the 
theory at large N c and at finite T and \i. There are two parameters we need 
to fix in vacuum: the 3-momentum cutoff A and the 4-Fermion coupling G. 
With f n ~ 92MeV for QCD, we have the relationship 

fl = K f ^e(A-lp|) ^ m| _ 



with the constituent quark energy E p = y\p\ 2 + Mq. On the right hand side 
of this equation, we have dropped all terms which vanish in the limit A — > oo. 
Note that f% should be of order N c and G ~ 1/N C for large A^ c . We also have 

(15) 

where the dynamical quark mass Mq is obtained as the solution of the gap 
equation 

Mj{ln(2A/M Q )-i} = A>(l-^y . (16) 
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For Nf = 2 taking m = 5 MeV with constituent quark mass Mq = 320 MeV 
and f n = 92 MeV x ^3/N c one gets: A = 646 MeV and G = 10.2 GeV~ 2 x 
(3/AQ. 

At finite temperature and density, the thermodynamic potential becomes 

m ») = U+ (M ~ G m)2 - T - £ / ^Tr In S-\(2n + 1)T, p)/T . (17) 
The trace in the above expression can be done explicitly leading to 



n = -2N f T 



d 3 p 



(2tt) 3 
f d 3 p 



Trln{l + Le~ {E ^ l)/T ) + (L -»■ L\ \l 

U . 



^e(A 2 - \pf) + (M m)2 



Taking $ = 1 in Eq. (lT8i) reproduces a standard NJL potential 



n 



NJL 



-2N c N f T 
- 2N c N f 



d 3 p 
J2^f 
d 3 p 



ln(l + e-^-ri) + ln(l + e ~^ +tM) ) 
(M - m) 2 



2G 



(19) 



which quantifies the thermodynamics of interacting quarks in the deconfined 
phase. 



4 Solving the Fermion Sector of the Theory at large N c 



The fermionic contribution to the thermodynamic potential is obtained from 
Eq. (USD as 



SSI* 



-2N f T 



d 3 p 

J2^y 



Trln{l + Le-K E »-" ) } + (L -> L 



7' 



(20) 



Suppose we are at /i < Mq. In this case the exponential factors inside the 
logarithm are all less than one. Moreover, L is a unitary matrix, LL^ = 1, so 
that we can always expand the logarithm as a series in L (To see this, work 
in a diagonal representation where each element of L is a pure phase). Using 
Tr lm4 = In detA for a matrix A, one finds 



8 



Trln(l + Le-^-ri) 



In 



1 + N c $ e-W'-ri + F 2 {(L 2 ), (L)) e- 2l3{E ?- 



■m) 



+ ■ ■ ■ + F p ((L p ), (LP- 1 ), • ■ ■ , (L» e ~^(^) + ■ ■ ■ + e" 



N c /3(E p -n) 



:2D 



The coefficient of p-quark contribution F p is a function that contains the trace 
of at most L p . 

In the confined phase, the first nonzero contribution occurs when the de- 
terminant is expanded to order L N °. This contribution is of order e~ KNc for 
temperatures of order Aqcd and fi a finite amount below M. Therefore quarks 
do not contribute to the effective potential in the confined phase, since they 
are exponentially suppressed. In the de-confined phase, all terms of order L p 
contribute, so the fermion are important again. Note however, that to compute 
the contribution of the fermion determinant requires evaluating contributions 
of order (L p ) which cannot be simply re-expressed in terms of (L). In fact 
to really deal with the large N c limit requires an effective potential for the 
Wilson line in the pure gauge sector which retains contributions such as (L p ). 
Nevertheless, even when in the de-confined phase, where in large N c there is an 
expectation value of L, it should be a good approximation to expand out the 
fermion determinant to first order in L. This is known for a free fermion the- 
ory where the Boltzman statistics result is accurate within 20% — 30%. Thus, 
to include the non-leading order effects in N c , we will therefore expand the 
determinant to first order in L. Note that in the confined phases of the theory, 
this should always be a very good approximation since there, in leading order 
in 1/N C , $ = and non-leading effects generate only a small expectation value 
for $. The terms of higher order in L are expected to be further suppressed. 

Consequently, when expanding the fermion determinant in positive powers of 
e -P\,E-n) £ or ^ < we conc i uc i e that whenever we are confined then baryons 
are exponentially suppressed. In large N c , there is no affect of temperature 
on the boundary between the quarkyonic and confined phase. There is no 
contribution of fermions to the expectation value of the a field, so that the 
chiral symmetry is unaffected by an increase in either density or temperature 
while in the confined phase. There is no feed back in large N c of the gluons onto 
the expectation values of the a field, again in the confined phase. We expect 
that the above is changed when going to finite N c , and that the expectation 
value of the a field will weaken as we go to higher temperatures. This is 
the conventional picture that increasing temperature destabilizes the chiral 
condensate. Note also that at finite N c , we should expect that the fermions will 
feed back upon the gluon potential. Assume N c is large so that we can expand 
the fermion contribution and keep only the first term in $, then this term will 
act to destabilize the first order deconfinement transition. Its effect increases 
as the baryon number density increases. We are therefore led to a picture like 
that of Philipsen and de Forcrand: The confinement transition weakens with 
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increasing baryon number density. Therefore, if there is a critical endpoint 
for realistic N c , it is the critical end point of the chiral phase transition, not 
that of confinement. We will discuss these non- leading effects in N c in a later 
section. 

Now let us determine the properties of the system when [i > Mq. We assume 
we are in the confined phase. We have to rearrange the determinant whenever 
fi> E p . We first note that the 

Tr In L = , (22) 
because it is an element of the SU(N C ) group so that we can rewrite 



Trln(l +Le~^ E ^) 

= TrlmLe-^-^ fl + lV^"^ 



= TrlnL + Trlne"^"^ + Trln (l + L^e~^~ Ep) ) 

= Tr pQi, - E p ) + Tr In (l + LV&^W) . (23) 

Thus, Eq. ( KOj) becomes 



6Q f = -2N f T J ^Tr [e(^ - fi) {ln(l + Le^-M)) + (L, /i - L\ -p) } 
+ 6(/i - E p ) {(3{v - E p ) + ln(l + tfe-^- E ^) + Qi -> -//)}] . (24) 

The second term in the above equation is part of the ideal gas contribution 
for a zero temperature degenerate gas of quarks. Note that it is no longer 
exponentially suppressed. Thus, non-interacting quarks contribute to the free 
energy. This contribution persists even in the confined phase. In the confined 
phase in large N c , this is the only term present, so it represents the contribution 
of quarkyonic matter. 

For large but finite N c if we are at low temperatures, far away from the decon- 
fined phase, the expectation value of $ is small. Therefore, we can expand the 
determinant to first order. This shows how we generate an explicit expectation 
value, and should show that its effects on the quarkyonic phase boundary are 
small. 

So lets assume that the expectation value of $ is very small in cold matter. 
Then, the vacuum contribution plus the finite chemical potential contribution 
of the fermions are of the form 
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-2N c N f J - 
-2N c N f v J 




d 3 p 



E p (Q(A-E p )-e(fi-E p )) 



quark 



d 3 p 



e(fi-E p ). 



(25) 



(2tt) 3 



The extremization of this effective potential leads to the gap equation, 



where the integral is cut off at the Fermi momentum pf = V f 12 ~ M 2 . The 
chiral phase transition point is determined by the vanishing of the chiral con- 
densate, M — > 0. Using Eq. (TIB"]) and the above gives us the solution for the 
position of the chiral phase transition as 



with the vacuum quark mass Mq. Putting in numbers, one finds that /i c hi ra i = 
1.04 M Q iniV c = 30. Note, that this equation tells us that the chiral transition 
occurs at densities slightly greater than that where the quarkyonic transition 
occurs (/j, = Mq). The cutoff A is independent of N c for large N c and thus the 
coefficient ln(2A/Mg) — 1/2 has only a weak iV c -dependence. The significance 
of this conclusion is subject to many uncertainties. One cannot therefore with 
certainty conclude nor rule out, that there is an intermediate phase of uncon- 
fined constituent quarks, as has been argued by Feinberg et al. [21]. In any 
case, the position of the chiral and the quarkyonic phase occur at numerically 
close, or perhaps identical, values of \i. 



5 Gluon effects on the chiral phase transition in large N c 

Let us begin by considering the response of the Fermion determinant to a non- 
zero expectation value of $. If we are at high temperature in the deconfined 
phase, then it is a good approximation to set $ = 1. In this case, we have an 
ideal gas of quarks whose thermodynamics is described by a pure NJL model 
given in (Tl9l) . In the following, we will be interested in the restoration of chiral 
symmetry where near the restoration point, we may approximate the mass 
as small. In this limit, one determines the chiral critical line from the gap 



The model with present parameters shows a first order phase transition in low 
temperatures. Nevertheless, an actual numerical calculation with m = 5 MeV gives 
almost identical critical value for a first order transition, // c hirai = 1.05 Mq. 




(26) 



A&taiCT = 0) = M% |ln(2A/M Q ) - i} , 



(27) 
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equation 



N c NfGA 2 GN c Nf 



2tt 2 



7T 



dpp 



1 1 

+ 



1 + ePte-v) i 4- e Pip+v) J ' 



(2? 



where the integration can be carried out analytically. Using Eq. (|T6|) for the 
constituent quark mass Mq in vacuum, this equation becomes 



f + yT 2 = Mj {ln(2A/M Q ) - \ 



(29) 



If we put in numbers, we see that everywhere in the de-confined phase, in the 
limit L = 1, chiral symmetry is restored for T c hi ra i ~ 177 MeV at /i — at 
iV c = 3. 

If we include the effect of nonzero by expanding the fermion determinant to 
first order in $ (an approximation which is good to 20%-30% even for $ = 1), 
one obtains 



n = u 



(M -m) 



2G 

5Q f = -2N c N f T 
+ 0(/i - E p ) 



2N f N c 



d 3 p 

W) 



E p Q{X-\p\) + 5£l 



3^P 



i ■ 



P{E P + I i) 



(2vr) 3 

/i — E p ) + $ (e-0O*-^> + e -0O*+*W) } 



(30) 



(31) 



where a difference between expectation values of $ and $ at finite \x is ne- 
glected. This leads to the gap equations for M and $, 



AVA 2 + M 2 - M 2 In 



'A+ VA 2 + M^ 
M 



m \ 2tt 2 



4$MTcosh(/i/T)if 1 (M/T) + - M) 
-4$ / M dE\jE 2 - M 2 cosh((£ - //) /T) 

JM 



PfH-M In ^ ^ 



T 3 C 



'iV 2 



-6 2 (T) + 6 4 $ 2 + &e$' 



2N C N 



f 



cosh(/i /T)M 2 TK 2 {M/T) 

+ 6(// - M) dEVW^M 2 sinh((£ - //) /T) 

Jm 



(32) 



(33) 



One can see that the gap equations take a simple form assuming a second 
order chiral transition, 
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Fig. 5. The phase diagram for N c = 3 obtained in our model with the current quark 
mass m = 5 MeV. The solid lines indicate a first order phase transition while the 
dashed lines cross over transitions. The critical end point (CEP) is indicated by a 
dot on the chiral phase boundary lines. 



fi 2 + 4$T 2 
/i 2 + 2T 



ln(2A/M Q ) 



16AUV 



-T 



f 



(34) 



(35) 



Thus, in the large N c limit they are two gap equations which describe the 
quark and gluon sectors without any interference. Eq. (IMj) coincides with 
Eq. ( |27l) and thus the chiral phase transition is indicated by a straight line 
/Uchirait^) = /^chirai^ = 0). This also dictates the order of phase transition for 
any T . Eq. (135]) determines the first-order deconfmement transition described 
by $(T) and thermodynamics now depends only on T. 

Fig. [5] shows the model phase diagram for N c = 3. The chiral and deconfme- 
ment cross over lines are identified as a maximum of derivatives DM/ dT and 
<9$ / <9T, respectively. The chiral and deconfmement lines are almost on top in 
a wide range of /i. Near the critical end point (CEP) the chiral and quarkyonic 
transitions are strongly coupled and the CEP appears near the intersection 
of those boundary lines. The chiral phase boundary is influenced the decon- 
fmement transition, and may weaken the chiral transition and result in the 
appearance of a CEP. 

The evolution of the phase boundaries with N c is shown in Fig. [61 The chiral 
transition lines move to larger /x and approach the vertical line. Correspond- 
ingly, the CEP is also shifted to the right with iV c and eventually disappears 
in the N c — > oo limit. The coincidence of the chiral and deconfmement transi- 
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Fig. 6. The deconfinement and chiral phase boundary for various N c . The horizontal 
line describes the deconfinement phase boundary in iV c — ► oo. The vertical line 
indicates the chiral phase boundary in N c — ► oo. The solid line indicate a first order 
phase transition while the dashed lines cross over transitions. The symbols represent 
the chiral CEP for corresponding N c . 
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Fig. 7. The deconfinement transition lines for various iV c The horizontal line de- 
scribes the deconfinement phase boundary in ff c ^ oo. The solid lines indicate a 
first order phase transition while the dashed lines cross over transitions. The symbols 
represent the CEP asociated with Z(N C ) symmetry. 

tions is unaffected by N c . Both lines are equally shifted upward approaching 
the horizontal line To = 270 MeV characterizing deconfinement transition 
temperature in pure gauge sector. 

Increasing N c further, the cross over of deconfinement turns into a first order 
transition and a CEP associated with the Z(N C ) symmetry appears at finite 
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Fig. 8. The phase diagram of our model for N c = 3 and for N c — > oo. The horizontal 
and vertical solid lines indicate the deconfinement and chiral phase boundaries in 
N c — > oo. The vertical broken line indicates the second order quarkyonic phase 
transition in iV c — > oo. 

fi. This behavior is indicated in Fig. [7J This CEP appears as a result of quark 
interactions which makes the transition weaken. Thus the CEP disappears 
again in the large N c limit since quarks do not affect deconfinement. 

Fig. [8] shows the phase diagram for N c = 3 and for N c = oo. The model 
describes three distinct phases in the large N c limit which agrees with the 
argument made in [23] . (There may or may not be a fourth phase in a narrow 
range of chemical potential where there is baryon number but chiral symmetry 
is not restored. However, such a phase might also be an artifact of the model 
used in our calculation.) For realistic N c = 3 the order of phase transitions is 
changed due to the quark-gluon interference. Nevertheless, one sees a remnant 
of the phase structure in large N c along with a deformation of the boundaries 
including finite N c effects. 



6 Summary and Conclusions 

In this paper, we have shown, that the first order phase transition and asso- 
ciated critical end point are driven by the chiral rather then deconfinement 
transition. The lines of both these transitions sit nearly atop the cross over 
associated with the remnants of the first order quarkyonic phase transitions 
of the large N c limit. The critical end point itself may appear where the cross 
over from the quarkyonic and the deconfinement phase transitions intersect. 

It is interesting that within this theory, chiral symmetry breaking may occur at 
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temperatures below that of deconfinement. This has been addressed in several 
papers, and might be associated with parity doubling [33ll34f 35| . 

Clearly, chiral symmetry restoration in QCD may be more complicated than 
that which appears in PNJL type models. The Fermi surface of the quarky- 
onic phase is presumably associated with confined particles, and it may be, 
that there is a chiral condensate associated with this Fermi surface. Therefore, 
although chiral symmetry may be approximately restored at or very near the 
quarkyonic transition, the full restoration of chiral symmetry might require 
much high energy density, and might ultimately be associated with the decon- 
finement transition. 

We have also seen that the cross overs associated with deconfinement and the 
quarkyonic phase transition are remnants of the first order phase transition 
which occurred in the limit of large N c . Thus, also at finite N c , the lines of 
cross overs reflect the phase structure seen in the large N c limit. It is inter- 
esting, that the chiral critical end point typically appears at the juncture of 
the confinement and quarkyonic cross overs. This property may be more gen- 
eral and may appear beyond the PNJL model analysis. This is because, the 
large change of baryon number density associated with the quarkyonic phase 
transition may drive a first order chiral phase transition until the increase in 
the number of degrees of freedom associated with deconfinement destabilizes 
it, resulting in a critical end point. 

In summary, it is fair to say that the conclusions from the PNJL model analysis 
are at best suggestive of what the true structure of high baryon number density 
matter might be. 
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